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Accurate knowledge of the gravitational-wave (GW) signal from inspiraling compact binaries is essential to 
detect these signatures in the data from GW interferometers. With recent advances in post-Newtonian (PN) 
theory and numerical relativity (NR) it has become possible to construct inspiral-merger-ringdown waveforms 
by combining both descriptions into one complete hybrid signal. While addressing the reliability of such wave- 
forms in different points of the physical parameter space, previous studies have identified the PN contribution 
as the dominant source of error, which can be reduced by incorporating longer NR simulations. In this paper 
we overcome the two outstanding issues that make it difficult to determine the minimum simulation length nec- 
essary to produce suitably accurate hybrids for GW astronomy applications: (1) the relevant criteria for a GW 
search is the mismatch between the true waveform and a set of model waveforms, optimized over all waveforms 
in the model, but for discrete hybrids this optimization was not yet possible. (2) these calculations typically 
require that numerical waveforms already exist, while we develop an algorithm to estimate hybrid mismatch er- 
rors without numerical data, which enables us to estimate the necessary NR waveform length before performing 
the simulation. Our procedure relies on combining supposedly equivalent PN models at highest available order 
with common data in the NR regime, and their difference serves as a measure of the uncertainty assumed in each 
waveform. Contrary to some earlier studies, we estimate that ~10 NR orbits before merger should allow for the 
construction of waveform families that are accurate enough for detection in a broad range of parameters, only 
excluding highly spinning, unequal-mass systems. Nonspinning systems, even with high mass-ratio {q > 20) 
are well modeled for astrophysically reasonable component masses. In addition, the parameter bias is only of 
the order of 1% for total mass and symmetric mass-ratio and less than 0.1 for the dimensionless spin magnitude. 
We take the view that similar NR waveform lengths will remain the state of the art in the advanced detector era, 
and begin to assess the limits of the science that can be done with them. 

PACS numbers: 04.30.Db, 04.25.Dg, 04.25.Nx, 04.30.Tv 
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I. INTRODUCTION 

A network of gravitational-wave (GW) detectors is prepar- 
ing to achieve a remarkable scientific goal: the first direct de- 
tection of GWs. This will not only test the predictions from 
Einstein's general theory of relativity, it will also open a new 
window to the universe, revealing details of the population, 
composition and formation history of various astrophysical 
objects [IJ. One particularly interesting and promising source 
of detectable GWs is the inspiral, merger and ringdown of 
compact objects, such as black holes or neutron stars. 

An important contribution to the effort of detecting the sig- 
nature of coalescing compact binaries in the noise-dominated 
spectrum of a GW interferometer is the accurate modeling 
of the expected signals. Only with an entire family of these 
theoretically predicted template signals is it possible to filter 
large amounts of data taken from the interferometers. In a 
"matched-filter" search (see e.g. [21), these data are convolved 
with the model signals and if the agreement exceeds some 
predefined threshold one claims detection and further exploits 
theoretical predictions to estimate physical parameters of the 
binary system, such as component masses and spins. 

In the case of a binary black hole (BBH) with comparable 
masses, at least two different approaches are needed to de- 
scribe the full motion and radiated GW content from the sys- 
tem. Post-Newtonian (PN) theory is an asymptotic weak field 
approximation that treats black holes as point particles with 
a relative velocity v that is small with respect to the speed of 



light c (for details, see e.g., |J3] and reference therein). The 
standard PN formulation is based on expanding the relevant 
quantities (such as energy and GW flux) in terms of the small 
parameter v/c. Depending on the details of the expansion, re- 
summation and integration of the resulting differential equa- 
tions, different waveform models for the early inspiral are 
known, commonly denoted by TaylorTn (with n = 1 , . . . , 4) 
OHEl, TayIorF2 |9-12| and TaylorEt |13"l45. A further in- 
spiral waveform family is obtained by mapping the two body 
problem to an effective one body (BOB) system with the ap- 
propriate potential |l5l, fTsLTTTl . 

All these analytical approximations break down in the 
strong gravity regime, and one has to perform computation- 
ally expensive numerical calculations in full general relativity 
to describe the complete dynamics. Since 2005 |18-20| sta- 
ble numerical-relativity (NR) simulations of a few orbits plus 
merger and ringdown to a final Kerr black hole have become 
a standard tool to consistently predict the last stages of a BBH 
coalescence ||2TI|. The exploration of the whole parameter 
space, however, has just begun and, for instance, long simula- 
tions of systems with mass-ratios higher than q — m2/m\ ~ 5 
are still exceptionally time-consuming. For current overviews 
of the field see 1 22 - 25 1 . 

An obvious goal is to combine PN and NR results to pro- 
duce "complete" waveform models. Such signals contain 
physical information up to frequencies higher than the pure 
PN templates, which becomes increasingly impoitant when 
the total mass of the system increases. According to re- 



cent studies, binary neutron stars as well as mixed black 
hole/neutron star binaries can be detected well by point- 
particle PN templates 1261 l27l assuming the current and antic- 
ipated performance of the Laser Interferometer Gravitational- 
wave Observatory (LIGO) 1 28 J . We therefore focus on com- 
plete waveform models for BBH coalescences in this paper 
By including systems with small total masses in our analysis, 
however, we effectively consider the detection problem for a 
broad range of possible compact binary systems, although the 
extraction of all physical effects requires further modeling in 
the neutron star case. 

Several approaches have already been suggested to analyt- 
ically build template families that include all the stages that 
the BBH undergoes. The BOB family has been refined by 
adding extra parameters that cannot be determined by PN cal- 
culations but are fixed by calibrating them to highly accu- 
rate NR data. This combination of BOB and NR information 
yields a description of the entire coalescence process in the 
time domain |.29r.36J , often referred to as BOBNR. A differ- 
ent time-domain description based on standard PN expansions 
was presented in fyP\ as a step towards modeling generic spin 
configurations. In this paper, we consider the direct match- 
ing of standard PN waveform models to NR data. In this ap- 
proach, PN data are used up to some point in time or frequency 
and NR results are then taken to describe the remaining part 
of the waveform. A phenomenological fitting of these "hy- 
brid" waveforms can then be performed to obtain an analytical 
closed formula (in the frequency domain) which interpolates 
between the physical parameters of the hybrids I381 441L 

All these procedures are subject to ambiguities and errors 
that limit the applicability of the final waveforms. Here we 
will focus on the error due to the PN contributions to the hy- 
brids only. Previous work has shown that the uncertainties in 
the NR waveforms and in the hybridization procedure make 
a negligible contribution to the overall hybrid waveform error 
budget |41-43|, and as such we will estimate modeling errors 
on the basis of the dominating inspiral part of the waveform. 
In the absence of a well-defined notion of the PN error, how- 
ever, we have to account for it simply by considering two dif- 
ferent PN descriptions of the inspiral signal, which are equiv- 
alent to all known orders of their Taylor expansion. We then 
quantify the effect of this ambiguous part of the waveform by 
calculating an appropriately defined inner product ("match") 
between both choices. This data analysis-motivated measure 
leads directly to conclusions about how useful hybrid wave- 
forms are in the presence of an ambiguous PN part, or con- 
versely, what requirements have to be posed in order to model 
waveforms accurately enough. 

As an important application of this procedure we started 
in a previous paper |42| addressing the question of how long 
numerical waveforms have to be in order to fulfill the accu- 
racy requirements for a PN/NR hybridization. Our analysis 
of nonspinning binaries with mass-ratio q E [1,4] and equal- 
mass binaries with spins (anti-)aligned to the orbital angular 
momentum (with Xi = Si/m'2 < 0.5) lead to the conclusion 
that NR simulations of such systems should cover 5 to 10 or- 
bits to be used in hybrids that satisfy the minimal accuracy 
requirement for signal detection. Bor larger mass-ratios and 



larger spins, our results suggested that far longer numerical 
waveforms were required. 

However, that study was limited due to the following re- 
striction: The efficacy of a model in a search is determined 
by the best match between the true waveform and any wave- 
form in the search model. This best match (called the "fitting 
factor") should be calculated not only by comparing two can- 
didates but by maximizing the match over all of the physical 
parameters of the model. With access to hybrids from discrete 
points in the parameter space, we were only able to maxi- 
mize the match over the total mass of the binary, and so our 
results were a (possibly very) conservative estimate of wave- 
form length requirements. 

Bven stronger requirements were presented in a number of 
other studies, where no maximization was performed at all 
(except over the initial phase and time-of-arrival of the sig- 
nal), with the intention of determining the waveform length 
requirements not just for detection, but also for parameter esti- 
mation. With these more stringent requirements, MacDonald 
et al. B3l as well as Boyle ll44ll . concluded that NR wave- 
forms generally have to be much longer than currently possi- 
ble to produce hybrids sufficiently accurate for both detection 
and parameter estimation. In addition, Damour et al. B3]| 
presented a detailed comparison of phenomenological wave- 
form models |40 41 1 and a recent member of the BOBNR 
family |33|. As part of their approach they find that in partic- 
ular systems with higher mass-ratio {q> 10) can be combined 
accurately with a standard PN approximant in the frequency 
domain only if the NR waveform contains hundreds of orbits. 

In this paper, we study hybrid accuracy and NR waveform 
length requirements in the context of fully optimized mis- 
matches, i.e., fitting factors. Put differently, instead of quanti- 
fying the reliability of a single waveform with fixed physical 
parameters we ask how accurate the induced waveform fam- 
ily is at that point in the parameter space. To do this, we first 
simplify the nonoptimized match calculation, showing that it 
can be performed without full numerical waveforms. We then 
generalize our procedure to optimize the match with respect 
to physical parameters, and to then calculate the fitting factor 
that is necessary to make estimates of NR waveform length re- 
quirements that are meaningful for GW searches. By looking 
at the parameter bias between the best-match waveform and 
the target signal, we also gain some insight into the parame- 
ter estimation errors due to the uncertainties in the waveform 
modeling process. 

In the following sections, we will develop this procedure 
step by step, starting with the stringent assumptions previous 
results were based on and subsequently relaxing them until 
we reach the final result. After providing a mathematical def- 
inition of the (mis)match as our notion of error (Sec. |ll|, we 



show in Sec. Ill that the mismatch between two hybrids is 
determined by the PN uncertainty and the relative power be- 
tween the NR and PN parts of the signals. Thus, our accu- 
racy estimate requires only amplitude information in the NR 
regime and we derive how this can be incorporated, includ- 
ing the effect of possible time and phase shifts of the entire 
waveforms. (This is similar to the procedure developed by 
Boyle in I.44J . where instead of NR data, BOBNR signals are 



taken as "ersatz" waveform. We use a slightly more general 
approach by incorporating only amplitude information of the 
phenomenological model BTI .) Along the way we compare 
with previous results in the literature, showing that we fully 
agree on nonoptimized mismatch errors. 

When we finally optimize these mismatches with respect to 



physical parameters in Sec. IV we find that the corresponding 
errors for waveform families are much smaller than assumed 
so far. In particular, based on our estimates (that are mainly 
limited by the choice of PN families compared to each other) 
we conclude that NR simulations that cover ^10 orbits are 
probably acceptable for most astrophysical applications dur- 
ing the Advanced detector era. This includes nonspinning 
binaries (for which significant improvements in PN approx- 
imants are less likely in the next five years), where we ex- 
plicitly show that this relatively small number of NR orbits 
is sufficient up to at least ^ = 10, and with astrophysically 
reasonable restrictions even for ^ = 20 and above. We also 
adopt the view that, since typical simulations will be of com- 
parable lengths over the next five years, our focus here and 
in future work should not be on prescribing ideal (and unreal- 
istic) waveform lengths, but on determining the limits of the 
science that we can do with them. 



II. PRELIMINARY CONSIDERATIONS 

We shall address the question of accuracy of BBH hybrid 
waveforms in the following sense. Let 
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be the complex GW strain that combines the plus and cross 
polarization of the GW as the real and imaginary part, respec- 
tively. It is constructed from its PN description hp^ and the 
NR part /znr. We assume that the transition from /zpN to /znr 
is enforced at a single frequency 
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where h denotes the Fourier transform of h and /,„ is the 
matching frequency. Such a procedure can be employed in a 
direct Fourier-domain construction of the hybrid PTI . but it is 
also approximately true for time-domain hybrid constructions. 
In the latter case, the transition is carried out at a time t,„, 
where the instantaneous frequency is (o{t„,) — '-"^ — 2nfm. 
Then, for (|2]i to be true, we have to assume that 

1 . the transition frequency in the Fourier domain is equal 
to the instantaneous matching frequency calculated in 
the time domain; 

2. the signal at times t < t,„ only significantly affects the 
Fourier domain for / < /„, and t > t,„ correspondingly 
determines the wave for / > /,„. 

These assumptions are not trivial since the Fourier integral is a 
"global" transformation. However, it was shown that assum- 
ing such a stationarity is reasonable in a regime where both 



PN and NR are valid pTl and time- and frequency-domain 
construction methods lead to very similar results |,42l . 

The final hybrid waveform is subject to several errors, and 
we account for these errors here simply by the fact that one 
could have taken slightly different ingredients /ipN and /jnr 
for the same physical scenario. These could be different post- 
Newtonian approximants and numerical data from different 
codes or different resolutions. Denoting the different wave- 
form models by hi and /12, we calculate the mismatch 
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where 00 and to are relative phase and time shifts between the 
waveforms and \\h\\^ — {h,h). S„ is the noise spectral density 
of the assumed detector, * indicates the complex conjugation 
and [fx.fi) is a suitable integration range. is called the 
overlap (or match) of the two waveforms. Throughout this 
paper, we will follow the choices of our preceding work |42|, 
i.e., /i = 20Hz and S,, is given by the analytic fit of the de- 
sign sensitivity of Advanced LIGO ["391 . The upper integra- 
tion bound /2 is given by our waveform model, and we use 
/2 = 0.15/M, although the results do not depend sensitively 
on this value (M is the total mass of the binary). 

Broadly speaking, the mismatch indicates how "close" h\ 
and /i2 are. Smaller values for ^^ represent smaller errors in 
the waveform model, given that h\ and h2 are approximations 
of the same signal. Direct conclusions can be drawn from cal- 
culating the mismatch: If ^ is less than some threshold, we 
regard the final hybrid as accurate enough for the purpose in 
question. For a maximum loss of 10% of the signals in the 
detection process, we can accept a mismatch of w 3%, dis- 
regarding the addition from a discrete template spacing. If 
we account for the latter, one may decrease the accepted mis- 
match in the waveform modeling to 1.5% (see a similar dis- 
cussion in |42|) or even 0.5% as suggested in |46|. 

A generally more stringent requirement is that the uncer- 
tainty we have in the modeling is indistinguishable by the 
detector. Such a statement is obviously dependent on how 
"loud" the signal is in the detector. As discussed in |46| and 
further detailed in fl31l47l we can write the indistinguishabil- 
ity criterion as 



\\hi-h2\\^<e'^ 



(5) 



where the waveforms are optimally aligned in the sense of (Hh 
and e parametrizes the effective noise-increase due to model 
uncertainties. The minimal requirement for h\ and /12 to be 
indistinguishable is e = 1, although |45| argues that e ^ 1/2 
and probably less are more reasonable thresholds. Manipu- 
lating (J5]l under the assumption of equal norms leads to the 
equivalent inequality (see the calculation in ll48l ) 
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where Peff = |l/i|l/e is the effective signal-to-noise ratio (SNR) 
of the signal. 



When we later calculate ^ as a measure of the error in 
hybrid waveforms, we can set various thresholds based on 
^ < ^max or Eq. (J6]l to evaluate the reliability of current 
models. A potentially very useful application is then to con- 
clude which matching frequency is needed (i.e., how long do 
the numerical waveforms have to be) to ensure the desired ac- 
curacy. 



III. HYBRID MISMATCHES 

Having introduced the mismatch between supposedly 
equivalent waveform models as our notion of error, we shall 
devote this section to simplifying the mismatch calculation of 
two hybrid waveforms, optimized only with respect to a rela- 
tive time and phase shift. As we have pointed out in the intro- 
duction, this is not the complete procedure to assess the model 
accuracy in terms of signal detection because the optimization 
with respect to physical parameters is not considered. How- 
ever, we first need to develop some insights into this sim- 
pler procedure to eventually generalize it in the next section. 



overlap reads 



The results presented in Sec. Ill C therefore have mainly illus 



trative character, showing that our simplified approach fully 
agrees with previously published results, but it leads to overly 
conservative requirements, e.g., for the length of NR wave- 
forms. 



A. General procedure 

Before we calculate mismatches for many different scenar- 
ios, we establish a few more assumptions to gain some in- 
sights on the structure of Eq. (HI. These will allow us to pro- 
pose an approximation to the mismatch between two hybrid 
waveforms that can be calculated without the need for any NR 
data. 

In addition to (|2]), we further assume: 

3. following ll2ni4n we regard the error on the NR side as 
small, negligible compared to the uncertainties PN in- 
troduces up to currently practical matching frequencies. 

4. independent of the PN approximant that is used, the 
norm of the waveforms are to high accuracy the same 
(i.e., only the phase is affected). This is reasonable to 
take as a good approximation, because the amplitude 
description in PN is usually formulated as a function of 
the orbital frequency |49-51 1 (which we again identify 
with the content on the Fourier side as well) and the 
mismatch is much more sensitive to phase differences 
than to amplitude discrepancies. 

Let us now consider a BBH system with fixed physical pa- 
rameters. Our error measurement assumes the construction of 
two hybrid waveforms that differ in the PN part only. Their 
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where A, = |/i, | and 0, = arg/z,. The effect of a time and phase 
shift of one waveform with respect to the other is explicitly 
written out in the second exponential term. 

Assuming two PN models (PNl and PN2) combined with 
the same NR waveform we trivially obtain the phase differ- 
ence 



-<j>2 = 



0PN1 — ^N2 
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Note that (IHll is only true for one particular alignment of the 
two waveforms, any other relative shift in time or phase leads 
to an additional dephasing, also beyond /„,. Since we have 
separated this effect explicitly in (JTl), we are, however, free to 
write 01 — 02 as in dSll. The open question is the functional 
form of the PN phase-difference (or simply the PN phase er- 
ror) in the case where the NR part of h\ and h2 are perfectly 
aligned. Here we have to apply an actual matching proce- 
dure, although we can use any preferred method without hav- 
ing NR data at hand. The key property of (J8]l we are exploiting 
is that only PN-PN differences are taken into account, and a 
direct PN-NR comparison is not necessary. The only input 
we need from NR simulations is the amplitude \h\ —Ai = A2 
for / > /,„. A good estimate for that can be taken from 
phenomenological models, such as 1401 or BTI . where the 
Fourier-domain amplitude is approximated by a closed-form 
analytic description. A similar approach was recently sug- 
gested by Boyle l:44J who realized that it is sufficient to com- 
bine PN approximants with ersatz NR data which he takes 
from the EOBNR model |[l5][l6l[30l|3l- We independently 
derive an algorithm here that is based on the same perceptions 
but highlights that no NR phase information at all is needed. 

The final global time and phase shift used in dTli to max- 
imize the overlap is simply a (phase shifted) inverse Fourier 
transform of the remaining integrand. Its maximal real part is 
obtained by choosing 00 (for any to) such that the generally 
complex number lies on the real axis. 

Based on that, our final algorithm for estimating hybrid 
mismatch errors caused by the uncertainty in the PN model 
is the following 

1. Calculate the two different PN waveforms expressing 
the uncertainty to be quantified. 

2. Apply the matching procedure such that one PN ap- 
proximant is matched to the other at /,„ (as if it were 
the NR waveform). 

3. Fourier transform the aligned PN waveforms and keep 
the data for /G [/i,/,„]. 

4. Complete the waveforms in the Fourier domain by us- 
ing an existing expression for the amplitude in the range 
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FIG. 1. Hybrid mismatches in the equal-mass, nonspinning case. 
Black solid lines are mismatches of actual TaylorTl/T4+NR hy- 
brids, whereas the dashed (red online) lines are our estimates ob- 
tained without directly using any NR data (NR amplitude taken from 
a phenomenological model). The gray dotted lines describe the PN 
mismatch contribution derived in Eq. ^ that does not include a pos- 
sible dephasing of the NR part. The matching frequencies for each 
set are from bottom to top Mco,,, = 0.04, 0.06, 0.08. 



/ > fm, e.g. from ll40ll4Tll or from a short NR simula- 
tion. Set the phase in this regime to (or any other 
function, but equal for both hi and /J2). 

5. Calculate the overlap of hi and /12 by maximizing the 
magnitude of the inverse Fourier transformation. 

To test the efficacy of our approach, we compare our esti- 
mate with the mismatch of actual hybrids consisting of either 
the TaylorTl or TaylorT4 approximant (in the form detailed 
in 1.41.1 ; see also Sec. [Ill Cl an d the numerical data from the 
SpEC equal-mass run 1521 I53l . The matching frequencies are 
chosen as Mco,„ = 2nMf,n € {0.04,0.06,0.08}, and the stitch- 
ing procedure is carried out in the Fourier domain as explained 
in Pn . The agreement illustrated in Fig. [T] is excellent in all 
cases. As expected by the relatively small effect of the am- 
plitude on the mismatch calculation, our method proves to be 
fairly robust with respect to the chosen amplitude description 
in the NR regime. In fact, the dashed lines in Fig. [Tluse the 
phenomenological model detailed by Santamaria et al. in BTI 
but there is no noticeable difference when we use the model 
presented by Ajith et al. in I.40J . 



B. Mismatch contributions 

The method presented above can readily be applied to es- 
timate the uncertainty of hybrids with the caveats mentioned 
at the beginning of Sec. Ill and we shall do so in Sec. |IIIC 
For now, however, let us manipulate the mismatch (Hh further 



to separate the various contributions to it. We make this im- 
portant aside to point out that, although only the PN contribu- 
tion is considered as ambiguous in our approach, its influence 
on the final waveform error is twofold: directly through the 
(power-weighted) PN mismatch and in terms of an additional 
dephasing, also of the "exact" high-frequency part. 



We can see these two effects separately through the follow- 
ing instructive lower bound on ^ which is obtained under the 
assumptions detailed above. 
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Here we introduced the notation \\h\\ 



i 1^ to specify the in- 
tegration range, ^pn is the mismatch of the PN part only, 
restricted to / < /„,. In the first line of (|9]l we use the fact 
that the amplitudes agree (in fact, we do not require point- 
wise agreement, only the norm is assumed to be the same) 
and that hi = h2 for / > /,„. The second line is a lower esti- 
mate because the maximization was originally carried out by 
shifting the entire waveforms relative to each other, whereas 
now we allow the maximization over the PN part alone. The 
final step involves the obvious relation \\h\\^ = \\h\\^r , ^ — 

The interpretation of (Ml is straightforward: the mismatch 
of hybrids is determined by the uncertainty of PN [restricted 
to the frequency range (/i,/,„)] multiplied by the fraction of 
power that is coming from the PN part of the wave signal. This 
fundamental error, independent of the actual PN/NR fitting, is 
directly inherited from the differences of standard PN approx- 
imants and any PN/NR matching cannot be better than the re- 
sult of Q. Therefore, one might think that analyzing the over- 
laps or fitting factors (or whatever strategy is appropriate) of 
different post-Newtonian approximants directly leads to con- 
clusions of how reliable the hybrid is for a particular choice of 
/,„. When we compare, however, the mismatch of actual hy- 
brid waveforms with the estimate (|9]) we find that the latter is 
considerably less than ./#. An illustration of that is included 
in Fig. [T] where we show the lower bound (J9]) in comparison 
with the actual (and accurately estimated) mismatches. 

Why is the hybrid disagreement that much greater than 
what is expected from PN in the given frequency range? The 
reason can be identified from the derivation of (|9]), where we 
effectively allow an optimal alignment (for each M) of both 
PN models while independently keeping the NR part per- 
fectly aligned. In a true hybrid mismatch calculation, one the 
other hand, a time and/or phase shift always affects the entire 
PN+NR hybrid, and an optimal alignment of one part leads to 
a dephasing of the other. This effect is not caused by an erro- 
neous matching, but an illustration of the fact that the optimal 
choice of to and (J)q in the sense of Eq. (J4|i is mass (frequency)- 
dependent for the PN models we consider 

Finally, by considering the obvious generalization of (|9]l. 
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we can identify the three main contributions to the hybrid un- 
certainty: The PN and NR error, each weighted by the power 



they contribute to the signal and the misalignment caused by 
the fact that in the hybridization procedure the PN wave is 
aligned at high frequency which is potentially different from 
the optimal alignment for lower frequencies. The procedure 



reduced, and we only use the mass-weighted total spin 

X^XimilM + X2m2/M (12) 



introduced in Sec. Ill A automatically takes the dominant PN 
error and possible misalignments (also of the NR part) into 
account. 



C. Application 

Now that we have established an algorithm to predict the 
full waveform mismatches, we can exploit the computation- 
ally cheap procedure and calculate ./# for many different 
physical scenarios. Our aim is to show how "reliable" the final 
combination of PN and NR waveforms is in different points 
of the parameter space, assuming that the physical parameters 
are fixed from the outset. 

First, let us highlight again that ideally, we are interested 
in the mismatch of the approximate waveform model to the 
true one. Since we cannot calculate the latter (which would 
also make the whole discussion pointless), we estimate the 
PN uncertainty by calculating the mismatch between different 
approximants. This can certainly be no more than a rough 
estimate since we are not aware of any principle that would 
guide us to which approximants at which PN order should be 
compared in order to obtain a well-defined notion of the PN 
error. 

To still reach some understanding of the uncertainty in cur- 
rently used high-order PN models we present the anticipated 
hybrid mismatches when approximants commonly denoted by 
TaylorTl, TaylorT4 and TaylorF2 are used. TaylorTl and T4 
are solutions of ordinary differential equations in the time do- 
main describing the adiabatic inspiral of the BBH on quasicir- 
cular orbits, whereas TaylorF2 is a frequency-domain descrip- 
tion based on the stationary phase approximation. Details on 
these approximants can be found, e.g., in |I7][8l and references 
therein. We mainly employ the equations presented in f4Tl, 
but with an updated 2PN spin-spin contribution from |51 1, see 
||54| for a collection of explicit expressions. Throughout this 
paper, we always employ the highest currently determined PN 
order, i.e., 3.5PN accurate phasing with spin contributions up 
to 2.5PN (and incomplete terms at higher order) and the 3PN 
amplitude expansion [,50 J including up to 2PN spinning cor- 
rections iBTI . 

As in the construction of phenomenological models, we 
restrict the parameter space to black holes with comparable 
masses and spins aligned or antialigned with the orbital angu- 
lar momentum of the binary L (with its unit vector denoted by 
L). Then, each spin can be parameterized by just one dimen- 
sionless quantity. 



Si L 

^2 



1,2, 



(11) 



where m, and Si are mass and spin of the individual black 
hole, respectively. By exploiting a degeneracy in the spins, 
as observed in 1551 l56l . the parameter space can be further 



and the symmetric mass-ratio 



n 



m\m2 



(13) 



to label the different physical setups. (In fact, in the following 
analyses, each point with fixed x is represented by Ji = X2 = 

X-) 

To assess how the accuracy of currently feasible hybrid 
waveforms varies in the parameter space, we apply the algo- 



rithm outlined in Sec. Ill A for different mass-ratios ranging 
from equal masses to 4:1, with spin magnitudes from —0.9 
to 0.9 in each case. For every pair {rjjX) one obtains mass- 
dependent mismatches in the form of Fig.fllthat generally in- 
crease with increasing matching frequency Mco,„. 

Several plots illustrating this behavior can akeady be found 
in the literature. Contour plots of the mismatch as a func- 
tion of mass and matching frequency are the main result of 
Boyle |44|, and we obtain similar results by continuously 
varying Mco,,,, e.g., in Fig.[T] Taking the maximum mismatch 
with respect to the total mass instead (i.e., only considering 
the peaks in Fig. fT), Fig. 4 by Damour, Nagar and Trias in 
1 451 shows the inaccuracy of TaylorF2 hybrids compared to 
EOBNR as a function of the matching frequency. Fig. 1 1 by 
MacDonald, Nissanke and Pfeiffer in |43| presents a similar 
study with Taylor approximants and actual NR data. Given 
some slightly different choices in our approaches (especially 
lower cutoff frequency and detector noise curve) the results 
we obtain are fully consistent with the numbers presented in 
the articles mentioned. 

Generally, the conclusions ||43H45]| draw are sobering re- 
garding GW detections and parameter estimation. The mis- 
matches found are too high, current numerical relativity wave- 
forms are by far too short and hybrids are consequently too 
inaccurate. In the following, we illustrate the basis of these 
statements and expand the existing knowledge by exploring 
the parameter space. To reduce the dimensionality of the prob- 
lem, we calculate the maximum of the mismatch with respect 
to the total mass and fix McOm = 0.06 (which corresponds to 
10 GW cycles before the maximum of \h{t)\ in the equal-mass 
case). 

In Fig. l2]we show contour plots that compare either Tay- 
lorTl with TaylorT4 hybrids or TaylorTl with TaylorF2 hy- 
brids. The matching frequency is fixed at Mco,„ = 0.06. Cer- 
tainly, we could include many more variants of PN approxi- 
mants (including different versions of EOBNR), but we find 
it sufficient to present some general conclusions that become 
already clear from the examples chosen here. As reported be- 
fore [42, .451 we see that deviating from equal-mass cases, the 
disagreement generally becomes larger. This effect is even 
more pronounced when increasing spin magnitudes are con- 
sidered. Heuristically we can understand the worse perfor- 
mance for increasing spins by the simple fact that spin con- 
tributions are only included up to 2.5PN order, whereas non- 
spinning terms are known up to relative 3.5PN order. Surpris- 
ingly, the 'island' or 'band' of minimal mismatch does not oc- 
cur strictly around vanishing spin magnitudes, indicating that 
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FIG. 2. Contour plots of the mismatch (in %) between different fictitious hybrids, as a function of the symmetric mass-ratio r\ and equal 
aligned spins with dimensionless magnitude X- Left panel: PN part either defined by the TaylorTl or TaylorT4 approximant. Right panel: 
Comparison of TaylorTl and TaylorF2 in the PN part. 
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3.93 (23) 
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25 


49 
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0.75 (730) 
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TABLE I. Faithfulness of hybrid waveforms based on a Tay- 
lorTl /TaylorF2 comparison for selected physical parameters. The 
required matching frequency is reported if either a 3% maximal mis- 
match ./# can be tolerated or if the error should be indistinguishable 
for SNRs less than 20, see |6l. The parentheses indicate the number 
of GW cycles to the maximum of |/i(/)|. The two right columns as- 
sume MoJm = 0.06 and give the minimal mass, where the waveforms 
are accurate enough in the sense described above. 



different approximants can by chance agree extremely well in 
some portions of the parameter space. For completeness, let 
us report that the TaylorT4/TaylorF2 mismatch yields a pat- 
tern similar to the right panel of Fig.l2]but with minimal values 
moved to weakly positive spins. 

The conclusions suggested by Fig. l2] and results from pre- 
vious work 141, 43 , 44 1 are indeed disappointing. If the mis- 
matches caused by different PN approximants actually rep- 
resent a reasonable estimate for the uncertainty in currently 
practical hybrid waveforms, then values up to ..# w 50% are 
certainly unacceptable. Reducing the matching frequency, 
thereby demanding longer NR waveforms, does reduce the 
mismatch everywhere, but it leads to unrealistic requirements 
in many portions of the parameter space. 

To illustrate this, Tablellladdresses two important questions 
by analyzing the TaylorTl /TayIorF2 hybrid mismatches in se- 
lected points in the parameter space. First, what is the required 
matching frequency if a desired accuracy has to be fulfilled? 



Note that due to our algorithm we overcome the restriction of 
currently available NR waveform lengths that the authors in 
BTl |431 were facing. We also do not rely on assuming a par- 
ticularly promising "candidate waveform" to act as a long NR 
waveform as was done in 1 42 44 1. In fact, phase information 
above Mo,„ is not required and does not enter the result; we 
can simply apply our algorithm to arbitrarily small matching 
frequencies. For each set of parameters we maximize the mis- 
match with respect to the total mass M (which we, however, 
restrict to M > SM© for computational reasons) and thus ob- 
tain the monotonically increasing function nmxM ^ {M (Om) ■ 
By demanding either ^ < 3% as the most relaxed require- 
ment or the more stringent case of indistinguishable differ- 
ences for effective SNRs of at most 20 [see (|6|] we obtain 
the values given in Table 11] In parentheses we also give the 
number of gravitational-wave cycles from d(l)Gw/dt — cOm to 
the maximum of |/i(f)| as predicted by the phenomenological 
waveform model |41 1. 

It is unlikely that the typical length of "long" numerical 
waveforms will change by an order of magnitude before the 
advent of Advanced LIGO, and so a more practical question 
is: given a currently achievable NR waveform length, in which 
mass-range is the PNh-NR hybrid accurate enough? As an ex- 
ample we assume again a matching frequency of McOm = 0.06 
and show on the right-hand side of Table|l]the minimal masses 
the hybrid is accurate for in the sense detailed above. For 
comparison, the pure NR part occupies the entire frequency 
band down to 20Hz for masses M > 91 Mq. Note that, distinct 
from 1 45 1, we do not consider errors above Mco,,, since we are 
concerned with hybrids and not possibly fitted closed-form 
waveform models that introduce additional errors. Therefore, 
our values for Mmin are less than the corresponding results in 
145 1 that are based on the comparison of EOBNR and the phe- 
nomenological model of POl . 

The obvious message from Table 11] is that in general ex- 
tremely long NR simulations would be needed to overcome 
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the intrinsic uncertainty in standard PN formulations for given 
physical parameters. For NR waveforms containing so many 
cycles our assumption that their intrinsic error can be ne- 
glected is possibly no longer valid, which would lead to even 
higher modeling errors. Anyway, the numbers presented are 
only an "order of magnitude" estimate in this most conserva- 
tive approach. The reader should always keep in mind that 
our notion of error is based on comparing different, at high- 
est available order consistent PN descriptions and especially 
concrete statements for particular points in parameter space 
may be spoiled by an (un)fortunate choice of approximants 
(see a similar discussion in |42|). More importantly, as we 
shall show in the next section, fixing the physical parameters 
of the waveforms from the outset greatly overestimates the 
uncertainty for signal detection. 



IV. FITTING FACTORS 



A. General 



The accuracy assessment presented in Sec. Ill only allows 
for very limited conclusions about the actual utility of hybrid 
waveforms in various applications. Apart from the restric- 
tions coming from our limited understanding of the PN error 
there is also an important fact we have neglected so far: in as- 
trophysically relevant applications the knowledge of physical 
parameters like total mass, mass ratio and spin is never exact. 
If a set of hybrid waveforms constitutes a waveform family 
which is used to extract information from an unknown signal, 
then the standard matched-filter procedures rely on varying 
(and maximizing with respect to) such parameters. The accu- 
racy of the predicted "best-fit" parameters is once again lim- 
ited by the detector noise and the modeling error and even if 
the latter exceeds the first, one may still argue that a tolerated 
bias does not significantly reduce the scientific output from 
GW detections. 

In this section we shall therefore consider combinations of 
NR data with a particular PN approximant as the ingredients 
of an entire manifold of waveforms, parametrized by an abso- 
lute time and phase scale (fo and ^o) as well as the physical 
parameters introduced before: M (total mass), 7] [symmetric 



mass-ratio ( 13 1] and x [spin combination ( 12 1]. The efficiency 
of detecting a signal defined by to,(j)(j,M,rj and x is properly 
quantified through the fitting factor 



FF = max & 

M',ri',x' 



h,{M',n',x')M{M,n,x) 



(14) 



Note that the maximization with respect to fo and 0o is already 
included in the definition of the overlap (?, see (J7|. 

The accuracy threshold for detection we quoted before is 
indeed defined including this additional maximization, i.e., in 
terms of 



^FF = 1 — FF , 



(15) 



the detection of the target signal 1x2 fSl. The results in Sec. Ill 
are only a lower bound on this effectualness. 

The accuracy requirements for parameter estimation are 
naturally more demanding than those for detection. In the re- 
cent literature BTl l43l l45l l46l the faithfulness of waveforms 
was usually defined by the criterion (|5]l (without optimiza- 
tion with respect to physical parameters), thereby demanding 
that the maximal information can be extracted from the data 
without being restricted by the model itself Here, however, 
we want to understand faithfulness in the original sense in- 
troduced in |5| that is based on the difference of the target 
waveform para meter A with the recovered model parameter A 
for which (14 1 is maximal. If this bias AA ::= A — A is small 



enough, we can still accept the waveform model family as suf- 
ficiently accurate, even for parameter estimation. Therefore, 
by analyzing ^ff and the corresponding parameters we can 
sensibly make analogous conclusions as before, but based on 
the actual optimization strategy that is employed in current 
template-based GW searches. 

Because of the additional freedom of varying physical pa- 
rameters we now have to calculate the ambiguity function 



s^[X\X)^ff[h{X')M{^)\ 



(16) 



between hybrids constructed from the same set of NR wave- 
forms but members of different PN approximants. It depends 
on the parameters of the waveforms, X and X, as well as the 
waveform models themselves. 

Since the phase difference above Maim in the overlap in- 
tegral ([TJi does not vanish generally for A ^ A, we have to 
slightly modify the algorithm presented in Sec. [Ill A In partic- 



If a waveform family {hi] satisfies ^ff(/ii,/!2) < ^max 
(with sufficiently small ^max) then it is said to be effectual in 



ular, we now need an estimate of how small changes in phys- 
ical parameters affect the phase difference in the assumed NR 
regime. (The PN regime is affected as well, but there is no 
qualitative difference to the PN comparison incorporated be- 
fore.) One possible strategy to quantify phase changes along 
variable physical parameters is to perform a number of nu- 
merical simulations and interpolate between the data obtained. 
Depending on the density of samples in the Tj and x directions 
(the scaling with M is given trivially by a single simulation), 
such a procedure can be very time- and resource-consuming. 
However, the phenomenological fittings performed in 1*381 4411 
have utilized exactly this type of interpolation, and we conve- 
niently use the result of BTl here because the fitting there is 
localized to frequencies close to and in the NR regime. 

Finally, to ensure the proper relative alignment, our algo- 
rithm to calculate ^ for arbitrary (in practice small) varia- 
tions in all parameters is to match different PN approximants 
to a phenomenological waveform (phase and amplitude) that 
is used above Mo,,, resulting in a hybrid /j (/; M, Tj , ;i^ , fo , 0o ) • 

Let us highlight that although we are now building 
PN-nphenomenological hybrids our analysis is not assessing 
how accurate individual waveforms describe the entire coa- 
lescence process. Note for instance that we could have in- 
troduced this hybridization concept already in the previous 
section, but, as we have shown, the phase above the match- 
ing frequency did not enter the overlap calculation. Similarly 
now, we use the phenomenological phase description merely 
to model the M-, Tj- and ;if -dependence at higher frequencies. 
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FIG. 3. The ambiguity function \\i>\ between two phenomenological 
waveforms |41 1, where /12 is fixed with r\ = 0.2, % = 0.3 and the 
total mass as indicated in the plots. The parameters of h\ are varied 
individually while the others are kept constant at the values of hi, 
respectively. 



Figure I3] illustrates what kind of information we are using by 
plotting slices of the ambiguity function of the phenomeno- 
logical model with itself for the case r\ = 0.2 (mass-ratio 
w 2.6), X = 0.3 and M/M,,, e {10,50, 100}. In Sec. |m]we 
only exploited £/ —I for X — X whereas now we need an es- 
timate of the shape of £/ also for X ^ X (although for small 
\X'~X\). 



We can make two immediate observations from Fig. |3] Es- 
pecially for small masses we see that relatively small changes 
in, for instance, symmetric mass-ratio or total mass (the other 
parameters are kept constant, respectively) modify the wave- 
form considerably, so that the high mismatches for equal pa- 
rameters (reported, e.g., in Figs. [T] and l2]l could potentially be 
reduced drastically by only small variations in the physical pa- 
rameters of one model waveform. Although the formal crite- 
rion Q for faithfulness (or better indistinguishability) failed, 
the fitting factor could still be extremely close to unity with 
a minimal bias in the parameters. The second interesting ob- 
servation from Fig. [3] is that the width around the maximum 
of the ambiguity function increases towards higher masses so 
that a comparison of two waveforms is increasingly insensi- 
tive to parameter changes at higher frequencies. This in turn 
endorses our assumption that the fitting factors and biases we 
shall calculate are dominated by PN effects (and not the choice 
of data above Maim) for small masses, where the accuracy re- 
quirements turned out to be hardest to satisfy. 



B. Comparison with previous results 

Before exploring fitting factors across the parameter space, 
let us present two examples that illustrate the general con- 
clusions we shall draw in this paper We first come back 
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FIG. 4. The mass-optimized mismatch between equal-mass, non- 
spinning TaylorTl/TaylorT4+NR hybrids (black solid lines) com- 
pared to our NR-free estimate (dashed lines). The matching frequen- 
cies are McOm £ {0.04, 0.06, 0.08} from bottom to top. The gray lines 
show the results of nonoptimized mismatches for comparison, see 
also Fig. [T] The inset illustrates the relative bias in the total mass 
(matching frequencies in reverse order). 



to the canonical equal-mass, nonspinning case and the Tay- 
lorTl/TaylorT4 comparison that was employed before (see 
Fig.[T]and 1421 ). To test the validity of our approach we again 
compare our estimate to hybrids constructed with actual NR 
data (matched at Ma)„, G {0.04,0.06,0.08}, respectively). Be- 
cause of the unavailability of NR data with arbitrary 77 and 
%, we for now only maximize with respect to the total mass 
M. Note that the results shown in Fig. [4] fully agree with the 
analysis of Hannam et al. ||42l (see Fig. 6 therein). They not 
only confirm that our combination of PN and phenomenolog- 
ical data accurately predicts the disagreement of the "true" 
PNh-NR hybrids, one can also observe the striking improve- 
ment when the additional maximization with respect to M is 
taken into account. The peak mismatch without optimiza- 
tion was approximately 8.8%, 4.5% or 2.2%, depending on 
Mco,„. With mass optimization we instead find ^pp < 3.2%, 
2.0% and 1.5%, respectively. The relative bias in the total 
mass, (M — M)/M, is always less than 0.8% and the earlier 
the matching is performed the smaller the bias becomes. 

A subsequent question that has not been answered so far is 
to what extent further optimizations, say along the symmetric 
mass-ratio and the spin(s) of the model system, improve the 
agreement between the waveform families even more. Full fit- 
ting factor calculations are commonly used to compare wave- 
form models (see, e.g., 18] |45l), but they have not been em- 
ployed in the context of hybrid waveforms and studies of the 
required length of numerical waveforms. Reference |42| only 
applied a crude estimation of the effect an additional mass- 
ratio optimization has, and concluded that a (total) mass- 
optimization alone serves as a sufficient assessment of the full 
fitting factor We now find that this conclusion was incor- 
rect. We illustrate the effect of further optimizations through 
the comparison of TaylorTl- and TaylorF2 -based waveforms 
(matched atMo),,, = 0.06) in Fig. [5] The TaylorTl target sig- 
nal is fixed as a system with mass-ratio 4: 1 and spin % — 0-5, a 
point in parameter space that clearly fails all accuracy require- 
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FIG. 5. Mismatch between TaylorTl- and TaylorF2-based waveforms for mass ratio 4:1, x = 0.5 and matching frequency Mco„, = 0.06 (left 
panel). The mismatch is not optimized (top solid line), mass-optimized (dashed line) or optimized with respect to all physical parameters of 
the TaylorF2-based model (lowest dotted line). The bias in the parameters are provided on the right panel. 



ments when looking at Fig. l2] By maximizing with respect 
to M, however, the maximal mismatch drops from 32.2% to 
10.4%. Varying all three considered physical parameters fi- 
nally yields a curve with .^pp « 1.6% at maximum, making 
the TaylorF2-based family accurate enough for detection. The 
relative bias in the parameters are less than 1 % for M, of the 
order of 1% for Tj and < 10% for X- 

Note that a faithfulness analysis, as in Sec. Ill and Il43ll44l . 
would conclude that NR waveforms with many hundreds of 
cycles are necessary to produce hybrids (and consequently 
waveform models) that are sufficient for parameter estimation 
purposes. Here we see that waveforms that we might at first 
sight regard as far too inaccurate, in fact may yield relatively 
small parameter biases when embedded in a waveform family. 

The optimization algorithm with respect to physical param- 
eters is computationally more challenging than maximizing 
the inner product with respect to fo and 0o only. For each set of 
test parameters {ri,x) we have to construct a new waveform. 
Since TaylorF2 is an analytical closed-form PN description 
that is fast to evaluate and our matching to the phenomenolog- 
ical model is performed directly in Fourier space |41 1 we only 
consider TaylorF2 -hybrids as test waveforms hi. For the fixed 
target waveforms /12 we chose to employ the TaylorTl approx- 
imant, because it was shown in ll57ll that its (dis)agreement to 
premerger NR data is most robust over the considered param- 
eter space and |44J noted that a maximal uncertainty estimate 
involves comparing to TaylorTl-inspirals. 

Starting with equal parameters X = A, we search for the 
nearest local maximum of the overlap ^(/ii , h2) by varying A 
along the gradient of the overlap. Thus, we ensure a quickly 
converging improvement after a relatively small number of 
iterations. The results we present, however, do not take into 
account the entire distribution of the ambiguity function and 
are still only a lower bound on the fitting factor. Given the 
tremendous decrease in mismatch for relatively small changes 
in physical parameters we argue nevertheless that this local 
extremum should serve as a reasonable estimate of the error 
one has to assume in terms of the fitting factor. 

We repeated the exploration of the parameter space with a 
study similar to the one presented in Fig.l2] The matching fre- 
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FIG. 6. The maximum of the optimized mismatch (in %) for hybrids 
constructed either with TaylorTl (target signal) or TaylorF2 (tem- 
plate signal) and a matching frequency of McOm = 0.06. 



quency is again fixed at Mo),,, — 0.06 and we calculate ^pp, 
Eq. (fH), for masses 5Mq <M< 20Mq. We checked that 



the mismatch decreases towards the boundaries of this inter- 
val, so that the enclosed maximum can indeed be regarded as 
the global extremum. After performing this maximization of 
the mismatch with respect to M for fixed {rj;X)^ W2 present 
our results as a contour plot in Fig. l6] The structure is very 
similar to the pattern of the nonoptimized mismatch, cf. the 
right panel of Fig.|2] The obvious difference is, however, that 
calculating the detection-relevant quantity ^pF instead of the 
diagonal mismatch 1 — £/(A,A) results in numbers that are 
^10 times less than what was considered before as error esti- 
mates. 

This allows for very different conclusions: Even a moder- 
ate matching frequency like the one considered here leads to 
hybrids that are accurate enough for detection in a large por- 
tion of the parameter space. Simulating NR waveforms with 
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few (< 10) orbits should hence be good enough for many ap- 
plications considering systems with moderate spins and mass- 
ratios. Although this is a very broad statement, it is clearly 
distinct from previous analyses Il43lfl3l that concluded much 
longer NR waveforms are needed to sensibly connect them to 
standard PN approximants. 

Of course, Fig. |6] only shows the optimal agreement be- 
tween the two considered waveform families and one might 
fear that the difference between simulated and recovered pa- 
rameters is large in some parts of the parameter space. How- 
ever, as anticipated by Fig.|5] the bias in total mass and sym- 
metric mass-ratio are small, approximately ±1% and ±1.5% 
at most, respectively. The spin parameter x is uncertain by 
—0.15 < Ax < 0.05. A deeper analysis of these biases is be- 
yond the scope of this paper and results are likely more model- 
dependent than the general conclusions we present here. 

For completeness, we note that for increasing values of the 
simulated spin, Atj and Ax generally decrease from positive 
to negative values, AM increases at the same time. This corre- 
lation is expected from the form of the PN expansion, where 
modifications of M can be compensated at lowest order by 
changing rj inversely. Studies of PN approximants in 1 8 1 show 
similar tendencies, although the biases reported there are con- 
siderably higher due to the absence of a common NR part at 
high frequencies. The same holds for the comparison of com- 
plete models (including uncertainties in the NR regime) II45I . 
The modeling biases we find should be compared to statisti- 
cal errors of full waveform families. In the case of the non- 
spinning phenomenological model fS^l a Fisher matrix study 
as well as Monte-Carlo simulations were presented in |58|, 
and the uncertainties found for Advanced LIGO and signals 
of SNR 10 are AM/M < 3% and Arj/rj < 8% (M < lOOM©). 
These values are of the same order of magnitude as our results, 
and we take this as an indication that modeling errors do not 
vastly dominate the parameter estimation uncertainty. How- 
ever, further studies are underway |59| to determine statistical 
errors for spinning waveform models. 



C. Model accuracy for spinning systems 

These new results constitute much brighter prospects for 
currently feasible NR simulations than the conclusions drawn 



in Sec. Ill and I42 n45il . In certain parts of the parameter space, 
however, the mismatch error presented in Fig. l6] is still too 
high, particularly if one keeps in mind that gaining sensitiv- 
ity of GW detectors is extremely difficult on the hardware 
side and theoretical considerations should reduce this sensi- 
tivity as little as possible |45|. Therefore, .4^pp > 3% for 
highly spinning systems should be improved by considering 
lower matching frequencies. Equally important is the ques- 
tion of whether numerical simulations for systems with mod- 
erate spins and mass ratios can be considerably shorter than 
M(0,„ — 0.06 which we assumed so far. 

In Fig. It] we analyze the dependence of the mismatch er- 
ror by showing the maximum of ^ff as a function of Mco,„. 
We consider equal masses and mass-ratio 4: 1 with spins x G 
{0,0.2,0.4,0.6,0.8} in each case. Note that we do not include 
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FIG. 7. The fully optimized mismatch ./#pp as a function of the 
matching frequency Moim for equal-mass systems (solid lines) and 
mass-ratio 4: 1 (dashed lines). The considered spins in each case are 
;f e {0, 0.2, 0.4, 0.6, 0.8} from bottom to top. 



orbits 


equal-mass 


mass-ratio 4:1 


5 


1.5%: -0.76 <X< 0.60 
0.5%: -0.37 <;i;< 0.31 


3.0%: -0.95 <X< 0.55 
1.5%: -0.52 <x< 0.39 


10 


1.5%: - 1.00 <;t;< 0.70 
0.5%: -0.45 <X< 0.39 


3.0%: -1.00 <;;;< 0.68 
1.5%: -0.56 <X< 0.48 


20 


0.5%: -0.97 <x< 0.57 
0.2%: -0.28 <X< 0.22 


3.0%: -1.00 <;(;< 0.79 
1.5%: -0.92 <x< 0.55 



TABLE II. Range in spin parameter % where a given accuracy re- 
quirement (.^FF < 3%, 1.5%, 0.5% or 0.2%) is fulfilled. Each row 
specifies the assumed number of orbits before merger for the NR 
waveform (= number of GW cycles divided by 2). 



negative values of x here, because the fact that the mismatch 
error for ;if < is smaller and not monotonic in x (see Fig.l6]l is 
likely an artifact of our choice of PN approximants (recall the 
obvious differences in Fig.l2]i. As expected, Fig.lTlillustrates 
that reducing the matching frequency, e.g., from Mo,,, = 0.08 
to Mo,,, = 0.02, leads to an improvement in mismatch by a 
factor of 2 to 10, depending on the spin. 

Larger values of the spin generally yield larger mismatches 
which in turn leads to stronger requirements forMo,,,, assum- 
ing a given accuracy goal. This is unfortunate because the 
orbital hangup configuration of positive aligned spins deceler- 
ates the frequency evolution in the inspiral of the binary, de- 
manding even longer simulations for a given frequency range. 

As such extremely long NR waveforms may not be avail- 
able in the near future (including the Advanced LIGO era), 
we continue with a slightly different application of our re- 
sults: How reliable is a set of complete waveforms constructed 
with standard PN approximants and NR simulations covering 
5 (10, 20) orbits before merger (i.e., 10, 20 or 40 GW cycles 
prior to the maximum of |/i(f)|)? To quantify these uncertain- 
ties we have to combine an estimate of the minimal matching 
frequency allowed by such NR waveforms with the resulting 
mismatch error from Fig. IT] We calculate the first from the in- 
verse Fourier transform of the phenomenological model 1.4 fl 
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FIG. 8. The fully optimized mismatch of nonspinning target sig- 
nals employing either the TaylorTl or TaylorT4 approximant with 
model waveforms constructed with TaylorF2 inspirals. The assumed 
matching frequency is always Ma>,„ = 0.06. The bias in parameters 
is \AM\/M < 0.6% (0.16%), |At]|/77 < 1.0% (0.3%), \A%\ < 0.04 
(0.017), where the values in brackets indicate the restriction to g < 4. 



and the time derivative of the phase, Mw,,, w dsagh{t„)/dt, 
where arg(f„) — argh{t^ax) — n2n (n = 10,20,40, respec- 
tively), and fmax is the time of the maximum amplitude \h\. 
This spin- and T] -dependent value is then taken into the re- 
sults presented with Fig.l7]to estimate ^pp for each configu- 
ration. Note that we use a more pessimistic error estimate for 
antialigned spins {% < 0) by assuming the mismatches of\x\ 
due to the reasons discussed above. 

One kind of possible conclusion one can then draw is sum- 
marized in Tablelnlfor equal masses and mass-ratio 4: 1. Given 
an accuracy goal (which we take as either 3%, 1.5%, 0.5% or 
0.2%) we provide the range of spins in which hybrids with 
the specified number of NR orbits fulfill this goal. Note that 
the asymmetry in the spin parameter is only caused by the dif- 
ferent matching frequencies waveforms with constant length 
permit. Again, we can very clearly see that even relatively 
short waveform are good enough for detection. In fact, mis- 
matches of 0.5% are below the noise level for SNR 10, and 
differences of 0.2% are indistinguishable for SNR < 16 ac- 
cording to Eq. (|6]). However, one can also see from Table [ll| 
that doubling the number of orbits does not enlarge the accu- 
racy range dramatically in many cases, although such simula- 
tions would take far more computer power and time. 



D. Nonspinning unequal-mass systems 

So far, we refrained from explicitly calculating mismatches 
for mass-ratios > 4:1 here because our underlying phe- 
nomenological model was only calibrated to numerical sim- 
ulations with mass-ratios < 4: 1 . Pushing the model beyond 
these values would add another uncertainty in addition to the 
way we estimate PN errors akeady, and more elaborate studies 
(possible including different models such as BOl and variants 
of EOBNR) are needed to reach sound conclusions. 

Nevertheless, numerical simulations of higher mass ratios 
are potentially interesting, and we shall try to estimate their 
reliability on the basis of our (extrapolated) knowledge here. 



We restrict this study, however, to nonspinning target signals. 
These are the systems where we do not expect the PN errors 
to drop significantly on the timescale of Advanced LIGO (in 
contrast to spinning binaries, where higher-order PN terms 
may well be calculated in the next few years). We find that 
the agreement between TaylorTl- and TaylorF2 -based hybrids 
is exceptionally good along % — ^ (see Figs. l6] and [Sll. In 
contrast, the TaylorT4/TaylorF2 uncertainty increases towards 
higher mass-ratios (smaller values of 77) as we would expect 
from the form of the PN expansion. Therefore, we shall con- 
servatively base our statements on comparing TaylorT4 and 
F2 approximants in this section. 

To illustrate our argument, we plot in Fig. [8] the maximum 
of the fully optimized (i.e., with respect to M, 7] and %) mis- 
matches between TaylorF2 and either TaylorTl or TaylorT4 
hybrids, all matched to fictitious NR data atMco,„ = 0.06. The 
fixed target parameters are chosen as ;); = with the mass ra- 
tio q varying from 1 to 4 in steps of 0.5 as well as ^ = 10 and 
q = 20. While the comparison with TaylorTl yields weakly 
77 -dependent mismatches below 0.3%, TaylorT4 target signals 
exhibit a steeply increasing divergence from the model signals 
towards higher mass ratios. Its approximately exponential be- 
havior is well described by the following fitting formula 



logio^PP« -0.29- 14.1 77 



(17) 



which is included as a straight line in Fig. [8] A conservative 
estimate of the general model uncertainty would be the maxi- 
mum of both data series for each 7j, i.e., ( fTTj i for small rj and 
roughly constant ^pp w 0.12% for 7; > 0.1866 (q < 3). 

Evidently, a matching frequency of M©,,, = 0.06 is only 
good enough for 77 > 0.081 {q < 10.2) if a mismatch of at 
most 3% is tolerated. Again, reducing the matching frequency 
helps to increase the accuracy of the final waveform, and we 
systematically analyze how useful numerical simulations of 5, 
10 or 20 orbits before merger are in the nonspinning unequal- 
mass regime. For that, we calculate maxM^-#FF as a function 
of the matching frequency and the symmetric mass ratio, sim- 
ilar to what was done for Fig. IT] The matching frequency is 
then converted to orbits before merger as explained in the pre- 
vious section. 

In Table III we present our results in analogy to Table In] 
where we provided the range of the spin parameter x in which 
the waveform model meets certain accuracy requirements. 
Now we complement the picture by restricting ourselves to 
the nonspinning case; our error estimates are based on opti- 
mized TaylorT4/TaylorF2 hybrid mismatches, and we present 
the accuracy range in terms of the mass ratio. Note that, al- 
though only five orbits of NR data before merger are sufficient 
for detection for most of today's standard simulations (q < 6), 
even the computationally very challenging goal of 20 orbits 
before merger is not enough to reliably model mass ratios as 
high as 15 or more for arbitrary total masses of the binary. 

It should be pointed out, however, that we report the worst 
disagreement between the considered hybrids in the left col- 
umn of Table III i.e., we demand that the assumed accuracy 
requirement is satisfied for all values of the total mass. As 
discussed in |44| already, one should rather understand the 
mismatch error and the accuracy requirement as functions 
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orbits 


mass-ratio 


? 


= 20 


5 


3.0%: 
1.5%: 


q<%.9 
^<6.8 


maXM^FF^ 

2IM0 : 12%, 


a 15% (19Ms) 
63Mq : 0.3% 


10 


3.0%: 
1.5%: 


q< 11.4 
^<8.6 


max^ ^FF " 
2IM0 : 3.0%, 


:^8.2% (I3M0) 

63Ms : 1.6x 10"^ 


20 


3.0%: 
1.5%: 


q< 14.8 
q< 10.7 


maxM ^FF " 
2IM0 : 0.8%, 


:^5.7% (IIM0) 
63M0 : 6.4 x 10"^ 



TABLE ni. Accuracy of nonspinning hybrid waveforms, based 
on combining PN TaylorT4 or TaylorF2 data with NR waveforms 
of specified length (defined by the number of orbits before merge 
= number of GW cycles divided by 2). Left column: Range in 
mass-ratio where a given accuracy requirement (max^^FF < 3% 
or 1.5%) is fulfilled. Right column: Mismatch error for q = 20, both 
at maximum of all masses (location indicated in parentheses) and at 
astrophysically motivated minimal values of the total mass (see text). 



of the total mass. After all, binaries with larger total mass 
have higher SNR in the detector (for constant distance of the 
source). More important for us here is that some of the con- 
sidered astrophysical scenarios may not even exist or be ex- 
tremely unlikely, and if the modeling error exceeds accuracy 
thresholds in these regions, we do not have to bother 

We illustrate this argument with a concrete example: The 
(fictitious) waveform of a binary with mass-ratio 20: 1 exhibits 
the largest uncertainty at total masses less than 2OM0, de- 
pending on the matching frequency (the values for NR sim- 
ulations covering 5, 10 or 20 orbits before merger are given 
in parenthesis in the right column of Table [IIIji. If we only 
consider black holes as objects in the binary and follow obser- 
vational f60l and theoretical I16IJ evidence that their individ- 
ual masses are > SM©, then the lowest total mass to consider 
in our error analysis is instead 63M0. With our idealized as- 
sumptions, this is a regime where the mismatch drops mono- 
tonically with increasing total mass (due to the dominating 
amount of exact high-frequency data), and the maximal un- 
certainty at 63Mq proves to be more than sufficient for detec- 
tion purposes, even with only a few NR orbits; see Table III 
In this sense, modeling higher mass-ratios is more accurate 
than comparable masses, as |44] noted already for diagonal 
(nonoptimized) mismatches. 

One the other hand, one could argue that the smaller object 
in the 20: 1 -binary could also be a neutron star. If the compan- 
ion is a much heavier black hole, tidal effects are extremely 
weak [271 and the plunge is hardly affected from finite size 
effects of the neutron star ll62ll . Thus, we may hope to accu- 
rately capture these systems with a BBH template family as 
well, and smaller total masses have to be considered. Accord- 
ing to |63|, (proto)neutron stars are expected to have masses 
> \Mq, which is in agreement with current observations (see 
^A\ for an overview). Assuming the lower bound of IMq for 
the mass of a single compact object, we consequently have to 
consider total masses down to 2IM0 (for q = 20) which leads 
to higher modeling uncertainties in the waveform. However, 



yield a mismatch of only 0.8% at 21M0. Hence, even the the- 
oretically and numerically difficult unequal-mass regime may 
well be modeled with only a few NR orbits, given the astro- 
physical expected properties of such systems. 

Of course, these astrophysical limitations are highly uncer- 
tain, and the conservative error analyses are the ones presented 
in Table In] and the left column of Table III However, given 
that caveat, we conclude that currently feasible numerical sim- 
ulations are potentially good enough to model in combination 
with PN approximants an important fraction of the parameter 
space. 



V. DISCUSSION 

Predicting the GW signature of an inspiraling and merging 
BBH in General Relativity is inevitably associated with ana- 
lytical or numerical approximations to the full theory, which 
introduce errors in the final result h{t) or h{f). In this pa- 
per we estimated these errors by the distance between two ap- 
proximate solutions for each physical configuration. While 
neglecting uncertainties on the NR side, we assumed differ- 
ent standard PN approximants in a frequency range up to the 
point where the waveform is matched to an NR-based merger 
and ringdown model. 

We quantified the uncertainties by comparing the currently 
available 3.5PN (spinning contributions up to 2.5PN) versions 
of TaylorTl, TaylorT4 and TaylorF2 approximants. Introduc- 
ing a simple algorithm that only requires amplitude informa- 
tion beyond the matching frequency, we first confirmed pre- 
vious studies 1 42-44 1 that found that the mismatch error for 
fixed physical parameters greatly exceeds reasonable accuracy 
requirements, assuming typical NR waveform lengths. 

Instead of demanding extremely long numerical simula- 
tions to overcome this uncertainty in the modeling process, 
we refined the understanding of the waveform error by adopt- 
ing the actual data analysis strategy of detecting an unknown 
signal in noise-dominated interferometer data. In particular, 
assuming waveform families instead of individual waveforms 
naturally redefines the concept of distance by allowing physi- 
cal parameters to be varied in the mismatch calculation. 



as Table III shows, 10 NR orbits before merger would be vir- 
tually good enough for detection purposes, 20 orbits already 



The results presented in Sec. IV indicate then that the GW 
signatures for many astrophysically relevant systems can in 
fact be well modeled by straightforward combinations of stan- 
dard PN approximants and currently feasible NR simulations, 
covering < 10 orbits before merger. The accuracy has not 
yet reached a level such that detection and parameter esti- 
mation errors are limited only by the detector noise for high 
SNR events, and the intrinsic uncertainty of BBH models 
may exceed in some cases the anticipated deviations caused 
by non-black holes, making is impossible to identify them as 
such. Nevertheless, the reported disagreement among differ- 
ent BBH models and biases in the parameters are certainly 
tolerable for the first GW detections that are likely to have 
low SNRs (^ 10). While this is true for systems with mod- 
erate spins, one has to keep in mind that even our idealized 
setting yields mismatch errors for high values of spins that 
are of the order of a few percent, which increase for higher 
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mass ratios. Reducing the matching frequency poses unreahs- 
tic challenges for current NR codes, and either fundamentally 
different numerical approaches or advances in PN are needed 
to fully control the entire parameter space. 

While the next spin-contributions in PN theory may become 
available in the near future to further improve the modeling of 
spinning systems (see the recent calculations of higher-order 
spin-orbit contributions |65-67J), unequal-mass nonspinning 
contributions at 4PN order are unlikely to be calculated with 
established techniques soon. However, as we discussed for 
a binary with mass-ratio 20:1, astrophysical expectations are 
that such systems only form with a high total mass, thereby 
reducing the impact of PN uncertainties. Even for 20:1 bina- 
ries, our results suggest that NR simulations of less than 10 
orbits are sufficient. 

In summary, we found that not single hybrid waveforms, 
but rather the embedding in the waveform manifold, results 
in templates accurate enough for detection, even with today's 
limited number of NR orbits. The uncertainty in physical 
parameters we had to accept for this tremendous increase in 
overlap is rather small, ^ 1% in mass and symmetric mass- 
ratio, and ^0.1 at most for the spin parameter %. For nearly 
equal-mass systems, the individual masses of the constituents 
are then only reliable to 




10%, 



(18) 



and it has to be decided whether this is good enough for astro- 
physical studies. 

Of course, our results rely on a number of assumptions that 
are reasonable in the range where we apply them, but we shall 
collect and discuss their generalizations and limitations below. 

First of all, our analyses are meant to provide a general con- 
cept of how to deal with modeling errors, instead of giving 
final answers. Especially, as we emphasized throughout the 
paper, we do not address the question of how accurate a par- 
ticular waveform model is. The statements formulated here 
are based on selecting PN approximants that are compared 
with each other, and our choices were made to illustrate the 
order of magnitude one generally has to assume for our no- 
tion of error This can be taken as a conservative estimate 
for all currently exsisting combinations of analytical and nu- 
merical relativity, because even a remarkable agreement in the 
overlapping region of both approaches does not necessarily di- 
minish the uncertainty of many ambiguous choices that enter 
the modeling of (up to thousands of) GW cycles in the inspiral 
waveform. Nevertheless, one should keep in mind that a par- 
ticular PN (or EOB)h-NR combination can be much closer to 
the real waveform than estimated here, as well as the possibil- 
ity that the PN ambiguity at consistent 3.5PN order generally 
underestimates the true error in the signal description. 

Two further essential assumptions should be noted: We ne- 
glect both the error of the hybridization procedure and any un- 
certainties beyond the matching frequency. Both assumptions 
are well motivated by previous studies 1271 I4TH43JI . but care 
has to be taken when generalizing their validity. For instance, 
from Fig.l6]or Table [Il]one might be tempted to conclude that 



actually very short NR waveforms are enough for modeling 
equal-mass, hardly spinning systems. This is certainly true 
from our results if the matching to PN can be done unambigu- 
ously. However, if there are too few cycles to align PN and 
NR signals properly, different matching procedures may lead 
to very different results. This aspect was not treated here as it 
can be checked separately, and it should only affect the result- 
ing waveform for very short (< 5 orbits) NR simulations. 

The other key assumption, the presence of exact high- 
frequency data, implies another important aspect to our re- 
sults. Not only do we say that the error of the NR part of 
the wave is negligible (an assumption that could easily be 
dropped if the NR mismatch becomes significant) we also use 
waveform families that directly resemble PN/NR hybrids. In 
other words, the additional error that is introduced in the phe- 
nomenological fitting and interpolation process is not taken 
into account here. Again, this is an error that can be quantified 
separately, but it has to be taken into account when interpret- 
ing the comparison of different complete waveform models, as 
was done in P31 . We merely state the fact here that in princi- 
ple PNh-NR combinations constitute sufficiently accurate tar- 
get waveforms for the construction of template families. 

This work can be complemented in many different ways. 
One obvious, yet involved extension is the completion of the 
parameter space by allowing arbitrary spin orientations that 
cause additional precession dynamics. Some steps towards 
building such hybrids have taken place already 1*681 16911 . but 
a deeper understanding of the waveform structure has to be 
gained before an extensive error analysis like the present one 
can be performed. Similarly, this study was restricted to the 
dominant spherical harmonic mode as it is crucial to under- 
stand and quantify the errors here first. Nevertheless, a final 
waveform model would have to include higher modes as well, 
and the algorithm we presented should be easily adaptable to 
these cases. 

Implementing more PN approximants and repeating our 
analysis with pairwise comparisons of various flavors of PN 
and FOB will help to fully understand the spread of equiva- 
lent descriptions of the inspiral process. When more contribu- 
tions to PN expansions become available the present analysis 
has to be repeated, hopefully reflecting the enhanced knowl- 
edge of the analytical approximation. This is especially true 
for spinning binaries, where calculations of higher-order PN 
contributions are expected in the next few years. 

Also, work is already underway ||59l to extend previous 
work |58| and relate the parameter uncertainties found in this 
study to statistical errors that are inevitably present for sig- 
nals with a given SNR in the detector Only these results will 
allow for statements about how useful current waveform con- 
structions are for parameter estimation and if the uncertainty 
in recovered parameters is dominated by the detector noise or 
the waveform model itself 

Finally and most interestingly, one should address the ques- 
tion of what kind of physics can be achieved given a certain 
performance of complete waveform models and, of course, 
given real GW detections with the upcoming generation of 
interferometers. It will be particularly important to analyze 
whether a certain disagreement between signal and model can 
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be entirely explained by model uncertainties or if possibly un- 
known physical effects are the cause. This study serves as a 
first step to prepare for those kinds of questions. 
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